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After introducing the Szekeres and Lemaltre-Tolman cosmological models, the real-time cosmol- 
ogy program is briefly mentioned. Then, a few widespread misconceptions about the cosmological 
models are pointed out and corrected. Investigation of null geodesic equations in the Szekeres mod- 
els shows that observers in favourable positions would see galaxies drift across the sky at a rate of 
up to 10~ 6 arc seconds per year. Such a drift would be possible to measure using devices that are 
under construction; the required time of monitoring would be « 10 years. This effect is zero in the 
FLRW models, so it provides a measure of inhomogeneity of the Universe. In the Szekeres models, 
the condition for zero drift is zero shear. But in the shearfree normal models, the condition for zero 
drift is that, in the comoving coordinates, the time dependence of the metric completely factors out. 
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I. INHOMOGENEOUS MODELS IN 
ASTROPHYSICS 1 

Just as was the case with our earlier reviews UQ, 
we define inhomogeneous cosmological models as those 
exact solutions of Einstein's equations that contain at 
least a subclass of nonvacuum and nonstatic Friedmann 
- Lemaitre - Robertson - Walker (FLRW) solutions as 
a limit. The reason for this choice is that such FLRW 
models are generally considered to be a good first approx- 
imation to a description of our real Universe, so it makes 
sense to consider only those other models that have a 
chance to be a still better approximation. Models that 
do not include an FLRW limit would not easily fulfil this 
condition. 

This is a live topic, with new contributions appearing 
frequently, so any review intended to be complete would 
become obsolete rather soon. Ref. 1 is complete until 
1994, Ref. 2 is a selective update until 2009 and Ref. 
3 is a more selective update until the end of 2010. The 
criterion of the selection in the updates was the usefulness 
of the chosen papers for solving problems of observational 
cosmology. 

In the very brief overview given here we put empha- 
sis on pointing out and correcting the erroneous results 
that exist in the literature and are being taken as proven 
truths. Some of them have evolved to become research 
paradigms, with many followers, some others proceed in 
this direction. We hope to stop this process, which dis- 
turbs and slows down the recognition of the inhomoge- 



1 This article is based on the talk delivered at the Marcel Grossman 
Meeting 13 in Stockholm, 2012. 



neous models as useful devices for understanding the ob- 
served Universe. 

We first present the two classes of inhomogeneous mod- 
els that, so far, proved most fruitful in their astrophysical 
application: the Szekeres model J4J, and its spherically 
symmetric limit, the Lemaitre [5[ - Tolman [6[ (L-T) 
model. Then we briefly mention the real-time cosmol- 
ogy program, and we give an overview of the erroneous 
ideas. Finally, we present an effect newly calculated in a 
few classes of models that is possible to observe and can 
become a test of homogeneity of the Universe: the drift 
of light rays induced by nonsymmetric flow of the cosmic 
medium 



II. THE SZEKERES SOLUTION 

The (quasi-spherical) Szekeres solution [H, 0] is, in co- 
moving coordinates 



ds 2 = df 
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(2.1) 



where E(r), M(r), P(r), Q(r) and 5(r) are arbitrary 
functions and $>(t,r) obeys 



$, t 2 = 2E(r) 



3> 



(2.2) 



The source in the Einstein equations is dust, whose mass 
density in energy units is 
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Eq. (12.2[) implies that the bang time is in general 
position-dependent: 



J J2E + 2M/$ + ±A$ 2 



t-t B (r). (2.4) 



The general Szekeres metric has no symmetry. It con- 
tains the spherically symmetric Lemaitre [B[ - Tolman 
(L-T) model as the limit of (P, Q, S) being all constant. 
The latter is usually used with a different parametrisa- 
tion of the spheres of constant (t, r), namely 



ds 2 = dt 2 _ ^^dr 2 - R 2 (t, r) (dz9 2 + sin 2 M V 2 ) , 

' " lL (2-5) 
where R = $, still obeying (|2.2p . and (|2.3|) simplifies to 

2M, r 



R 2 R, r ' 
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The Friedmann limit follows when, in addition, 
r) = r5(t), 2E 1 = -fcr 2 where k = const is the FLRW 
curvature index, and is is constant. 



III. REAL-TIME COSMOLOGY 

There are ways in which the expansion of the Universe 
might be directly observed. The authors of Refs. [1, 
composed them into a paradigm termed real-time cos- 
mology. One of them is the redshift drift: the change 
of redshift with time for a fixed light source, induced by 
the expansion of the Universe. 

Consider, as an example, an L-T model with A, given 
by (|2.5j) with R obeying (|2.2[) . Along a single radial 
null geodesic, directed toward the observer, t = T(t , r) 
(where t is the instant of observation) the redshift is [Hj 



1 + z(t ,r) = exp 



R, tr (T(t ,r),r) 
VT+2W) 



dr 



(3.1) 



For any fixed source (i.e. constant r), eq. (|2.2j) defines 
a different expansion velocity R,t for A = and for A ^ 
0, and thus allows us to calculate the contribution of 
A to z via (|3.1j) . With the evolution type known, the 
expansion velocity depends on r, and so allows us to infer 
the distribution of mass along the past light cone. 

According to the authors of @, H[ , the "European Ex- 
tremely Large Telescope" (E-ELT) 2 could detect the red- 
shift drift during less than 10 years of monitoring a given 
light source. The Gaia observatory 3 could achieve this 
during about 30 years. 



2 Now in the planning, to be built in Chile. 

3 http:/ /sci.esa.int/science-e/www/area/index.cfm?fareaid=26 



For more on redshift drift see the contribution by P. 
Mishra, M.-N. Celerier and T. Singh in these Proceed- 
ings. 

The drift of light rays described in Sec. IVland following 
is another real-time cosmology effect. 



IV. ERRONEOUS IDEAS AND PARADIGMS 

The L-T model was noticed in the astrophysics com- 
munity - but: 

1. Many astrophysicists treat it as an enemy to kill 
rather than as a useful new device. (Citation from Ref. 
@: The Gaia or E-ELT projects could distinguish FLRW 
from L-T "possibly eliminating an exotic alternative ex- 
planation to dark energy"). 

2. Some astrophysicists practise a loose approach to 
mathematics. An extreme example is to take for granted 
every equation found in any paper, without attention be- 
ing paid to the assumptions under which it was derived. 

Papers written in such a style planted errors in the 
literature, which then came to be taken as established 
facts. In this section a few characteristic errors are pre- 
sented (marked by ^) together with their explanations 
(marked by ^ ) . 

^ The accelerating expansion of the Universe is an 
observationally established fact (many refs., the Nobel 
Committee among them). 

^ The established fact is the smaller than expected 
observed luminosity of the SNIa supernovae (but even 
this is obtained assuming that FLRW is the right cos- 
mological model). The accelerating expansion is an 
element of theoretical explanation of this obser- 
vation. When the SNIa observations are interpreted 
against the background of a suitably adjusted L-T model, 
they can be explained by matter inhomogeneities along 
the line of sight, with decelerating expansion fill. H3|. 

^ Positive sign of the redshift drift is a direct confir- 
mation of accelerated expansion of the space. 

'fc The sign of redshift drift is only related to acceler- 
ation when homogeneity is assumed[13j (see also Mishra, 
Celerier, and Singh in these Proceedings). 

^ Hq, supernovae and the cosmic microwave back- 
ground radiation (i.e. the size of the sound horizon and 
the location of the acoustic peaks) are sufficient to rule 
out inhomogeneous L-T models. 

These observations only depend on D(z) and p{z), 
and thus can be accommodated by the L-T model, which 
is specified by 2 arbitrary functions. Examples of such 
constructions are given in [H|,[l4j]. In fact, as follows from 
the Sachs equations, in the approximation of small null 
shear (which for most L-T models works quite well[15j), 
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there is a relation between D(z), p(z) and H{z), mean- 
ing that these 3 observables are not independent^^], and 
thus allow the L-T model to accommodate more data on 
the past null cone. 

^ The gravitational potential of a typical structure in 
the Universe is of the order of 10~ 5 and thus, by writing 
the metric in the conformal Newtonian gauge, one im- 
mediately shows that inhomogeneities can only introduce 
minute (« 1CU 5 ) deviations from the RW geometry. This 
also means that the evolution of the Universe must be 
Friedmannian (common argument among cosmologists). 

^ Even if the gravitational potential remains small, 
its spatial derivatives do not, and thus the model has 
completely different optical properties than the Fried- 
mann models. This was shown by rewriting one of the 
L-T Gpc-scale inhomogeneous models in the conformal 
Newtonian coordinates [la ]. The gravitational potential 
of this model remains small yet the distance-redshift re- 
lation deviates strongly from that for the background 
model. 

The question remains whether small-scale fluctuations 
(of the order of tens of Mpc) could also modify optical 
properties and evolution of the Universe. The problem 
is complicated as it requires solving the Einstein equa- 
tions and null geodesies for a general matter distribu- 
tion, which, with current technology, is not possible to 
do numerically. Therefore, the problem has been ad- 
dressed in a number of approximations and toy mod- 
els. Recent studies showed that the optical properties 
along a single line of sight can be significantly different 
than in the Friedmann model. Yet, if averaged over all 
directions, the average distance-redshift relation closely 
follows that of the model, which describes the evolution 
of the average density and expansion rate[l5|. Thus, the 
problem reduces to the following one: do the average 
density and expansion rate follow the evolution of the 
homogeneous model, i.e. is the evolution of the back- 
ground affected by small-scale inhomogeneities and does 
it deviate from the Friedmannian evolution? Some au- 
thors claim that matter inhomogeneities cannot affect 
the background and that the Universe must have Fried- 
mannian properties fl7l [Tsj . while others argue for strong 
deviation from the Friedmannian evolution fl9l |20| . Stud- 
ies of this problem within the exact models, like L-T, 
proved that under certain conditions the back-reaction 
can be large, while under others it remains quite small 
plj . leaving the problem unsolved. For an informative 
description of t he p roblem and techniques used to address 
it see Ref. (H,^. 

^ Fitting an L-T model to number counts or the 
Dl(z) relation results in predicting a huge void, several 
hundred Mpc in radius, around the centre (too many pa- 
pers to be cited, literature still growing). Measurements 
of the dipole component of the CMB radiation then im- 
ply that our Galaxy should be very close to the center of 
this void, which contradicts the "cosmological principle" . 



The implied huge void is a consequence 
of handpicked constraints imposed on the L—T 
model, for example constant ts- When the model is 
employed at full generality, the giant void is not implied 

a 

^ The cosmological principle is a postulate, 
not a law of Nature, it cannot say which model is 
"right" and which is "wrong" . 

^ The bang time function must be constant, oth- 
erwise the decaying mode of density perturbation is 
nonzero, which implies large inhomogeneities in the early 
universe (see, for example, Ref. [Hj|). 

'fc The bang time function describes the differences 
in the age between different regions of the Universe. In 
the L-T and Szekeres models it is also related to the 
amplitude of the decaying mode. However, these mod- 
els describe the evolution of dust and therefore cannot 
be extended to times before the recombination, when the 
Universe was in a turbulent state: rotation, plasma, pres- 
sure gradients all did affect the proper time of an observer 
(dr = dt goa(t, x 1 )). Eventually, even if the Universe 
started with a simultaneous big bang, by the time of re- 
combination, due to the standard physical processes, the 
age of the Universe would have been different at different 
spatial positions, giving rise to non-constant is(r) of the 
dust L-T or Szekeres model that takes over there. 

Moreover, the relation between the nonsimultaneous 
big bang and the decaying mode was established only 
for the L-T @, H3 and Szekeres [H( models. For more 
general models, not yet explicitly known as solutions of 
Einstein's equations, like the ones mentioned above, the 
connection may be more complicated and indirect. Thus, 
citing this relation for such a general situation is an il- 
legitimate stretching of a theorem beyond the domain of 
its assumptions (see also the next entry below). 

^ The L-T models used to explain away dark energy 
must have their bang-time function constant, or else they 
"can be ruled out on the basis of the expected cosmic 
microwave background spectral distortion" [27[ ■ 

*fc The papers that claim this parametrise their 
models with a set of very simple functions, which 
lack flexibility. This sets them on the wrong track from 
the beginning. 

^ In order to meaningfully test any cosmolog- 
ical model against observations, one must apply 
it at every step of analysis of the observational data. 
To do so, would require a re-analysis of a huge pool of 
data. C. Hellaby with coworkers [28] is working on such 
a program applied to the L-T model, but the work is far 
from being completed. 

Lacking any better chance, we currently use observa- 
tions interpreted in the FLRW framework to infer about 
the M{r) and ts(r) functions in the L-T model. This is 
justified as long as we intend to point out possibilities, 
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under the tacit assumption that these results will be ver- 
ified in the future within a complete revision of the ob- 
servational material on the basis of the L— T model. How- 
ever, putting "precise" bounds on the L-T model func- 
tions using the self-inconsistent mixture of FLRW/L-T 
data available today is a self-delusion. An example: the 
spatial distribution of galaxies and voids is inferred from 
the luminosity distance vs. redshift relation that applies 
only in the FLRW models. Without assuming the FLRW 
background, we know nothing about this distribution un- 
til we reconstruct it using the L-T model from the be- 
ginning. 

The L-T and Szekeres models cannot be treated as 
exact models of the Universe, to be taken literally in 
all their aspects. They are exact as solutions of Ein- 
stein's equations, but when applied in cosmology, they 
are merely the next step of approximation after FLRW. 
If the FLRW approximation is good for some purposes, 
then a more detailed model, when applied in a situation, 
in which its assumptions are fulfilled, can only be better. 

V. THE REDSHIFT EQUATIONS IN THE 
SZEKERES MODELS 

Consider two light rays, the second one following the 
first after a short time-interval r, both emitted by the 
same source and arriving at the same observer. The tra- 
jectory of the first ray is given by 

(t,x,y) = (T(r),X(r),Y(r)), (5.1) 

the corresponding equation for the second ray is 

(i, x, y) = (T(r) + T{r),X(r) + C(r), Y (r) + ^(r)). (5.2) 

This means that while the first ray intersects a hyper- 
surface r = r at (t,x,y) = (T,X,Y), the second ray 
intersects the same hypersurface not only later, but, in 
general, at a different comoving location. In gen- 
eral the two rays will intersect different sequences 
of intermediate matter worldlines. 

The same is true for nonradial rays in the L-T model. 

Consequently, the second ray is emitted in a different 
direction and is received from a different direction by the 
observer. Thus, a typical observer in a Szekeres space- 
time should see each light source slowly drift across the 
sky. How slowly will be estimated further on. As will be 
seen from the following, the absence of this drift is 
a property of exceptionally simple geometries (or 
exceptional directions in more general geometries). 

We assume that (CVO an d (d/dr)(r, (, ip) are small of 
the same order as r, so we neglect all terms nonlinear in 
any of them and terms involving their products. 

For any function f(t, r, x, y) the symbol A/ will denote 

f(t + t, r, x + C, y + VO - f(t, r, x, y) (5.3) 

linearized in (rX,ip). Note: the difference is taken 
at the same value of r. Applying A to the null 



geodesic equations parametrised by r we obtain the equa- 
tions of propagation of (r, (, ip) and (£, rf) = f (d/dr) ((, ip) 
along a null geodesic - see both sets of equations fully 
displayed in Ref. [Hj]. 

VI. REPEATABLE LIGHT PATHS 

There will be no drift when, for a given source-observer 
pair, each light ray will proceed through the same inter- 
mediate sequence of matter world lines. Rays having this 
property will be called repeatable light paths (RLP). 
For a RLP we have 

C = V> = £ = v = o (6.1) 

all along the ray. The equations of propagation of 
(r, C, ip, £, rf) become then overdetermined (3 equations 
to determine the propagation of r along a null geodesic), 
and imply limitations on the metric components. They 
can be used in 2 ways: 

1. As the condition (on the metric) for all null 
geodesies to be RLPs. 

2. As the conditions under which special null geodesies 
are RLPs in subcases of the Szekeres spacetime. 

In the first interpretation, the equations of propagation 
should be identities in the components of dx a /dr, and 
this happens when 

* d = $, fr -$,i$, r /$ = 0. (6.2) 

This means zero shear, i.e. the Friedmann limit. Thus, 
we have the following 
Corollary: 

The only spacetimes in the Szekeres family in 
which all null geodesies have repeatable paths are 
the Friedmann models. 

In fact, we proved something stronger. Since the drift 
vanishes in the Friedmann models, we have: 
Corollary 2: 

The presence of the drift would be an observa- 
tional evidence for the Universe to be inhomoge- 
neous on large scales. 

In the second interpretation, there are only 2 nontrivial 
(i.e. non-Friedmannian) cases: 

A. When the Szekeres spacetime is axially symmetric 
(P and Q are constant). In this case, the RLPs are those 
null geodesies that stay on the axis of symmetry in each 
3-space of constant t. 

B. When the Szekeres spacetime is spherically sym- 
metric (P, Q, S are all constant) - then it reduces to the 
L-T model. In this case, the radial null geodesies are the 
only RLPs that exist. A formal proof of this statement 
is highly complicated [29j . 

For non-radial rays in the L-T model, the non-RLP 
phenomenon was predicted, by a different method (and 
under the name of cosmic parallax), by Quercellini et 
al. [§1, Q . Their review contains a broad presentation of 
the real-time cosmology paradigm, oriented toward ob- 
servational possibilities. 
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FIG. 1: Density profiles and positions of the observers used in the numerical examples. 



VII. EXAMPLES OF NON-RLPS IN THE L-T 
MODEL 

The examples will show the non-RLP effect for non- 
radial null geodesies in two configurations of the L-T 
model, shown in Fig. HJ for different positions of the 
observers with respect to the center of symmetry. In Ex- 
ample 1 (see Fig. [2]) we use Profile 1, the observer and 
the light source at 3.5 Gpc from the center of the void, 
the directions to them at the angle 1.8 rad, and [30j 

ts = simultaneous Big Bang, 

p(to,r) — po [l + S — Sexp (— r 2 /cr 2 )] - the density 
profile at the current instant, 

r = R(to,r) - the radial coordinate, 

Po = p(t o ,0) = 0.3 x (3ff 2 )/(87rG) = 0,3p clitical - the 
present density at the center of the void, 

Hq = 72 km s _1 Mpc -1 - the present value of the 
Hubble parameter, 

S = 4.05, a = 2.96 Gpc. 

The source in Fig. [2] sends three light rays to the same 
observer, the first of which was received by the observer 
5 x 10 9 years ago, the second is being received right now, 
and the third one will be received 5 x 10 9 years in the 
future. Fig. [3] shows these rays projected on the space 
t = now along the flow lines of the L-T dust. The source 
is at the upper right corner of the graph and the observer 
is at the upper left corner. 

The time-averaged rate of change of the position of 
the source in the sky, seen by the observer is 

angle between the earlier and the later ray 
time interval = 5 x 10 9 years 




FIG. 2: The configuration of the light source and the observer 
in example 1. 

The rate of drift in the next figures is calculated in the 
same way. 

In Examples 2, 3 and 4 (Fig. HJ, the observer O is 
at Ro from the center; the angle between the direction 
toward the galaxy (^) and toward the origin is 7. For 
each 7 we calculated the rate of change 7 by eq. (|7.ip . 
and the graphs in Fig. [5] show 7 as a function of 7. All 
the examples have d = 1 Gly w 306.6 Mpc. Example 
2 (solid line) has Rq — 3 Gpc and Profile 1; Example 3 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
[rad] 



FIG. 3: Three light rays projected on the space of constant comoving time along the flow lines of the cosmic medium. Middle 
line: the ray received at the current instant; upper line: the ray received 5 x 10 9 years ago; lower line: the ray to be received 
5 x 10 9 years in the future. 



(dashed line) has Rq = 1 Gpc and Profile 1; Example 4 
(dotted line) has Rq = 1 Gpc and Profile 2 (for which 
5 = 10.0 instead of 6 = 4.05; this is a deeper void in a 
higher-density background). The amplitude is ~ 10~ 7 for 
(3) and ~ 10 -6 for (2) and (4). With the Gaia accuracy 
of 5 — 20 x 10~ 6 arcsec, we would need a few years to 
detect this effect. 



VIII. RLPS IN SHEARFREE NORMAL 
MODELS 

In the Szekeres models, the condition for all null 
geodesies to be RLPs was the vanishing of shear. This 
suggests that the cause of the non-RLP phenomenon 
might be shear in the cosmic flow. To test this supposi- 
tion, the existence of RLPs was investigated in those cos- 
mological models in which shear is zero 1311 - the shear- 
free normal models found by Barnes [32j . They obey 
the Einstein equations with a perfect fluid source and 
contain, as the acceleration-free limit, the whole FLRW 
family. 

There are four classes of them: the Petrov type D met- 
rics that are spherically, plane and hyperbolically sym- 
metric, and the conformally flat metric found earlier by 
Stephani [33[. In the Petrov type D case, the metric in 
comoving coordinates is 

ds 2 = (^\ 2 dt 2 -±(dx 2 +dy 2 + dz 2 ), (8.1) 




FIG. 4: The configuration for Examples 2, 3, 4. 



where F(t) is an arbitrary function, related to the expan- 
sion scalar 6 by 9 = 3/ F. The Einstein equations reduce 
to the single equation: 

w, uu /w 2 = /(«), (8.2) 

where /(it) is an arbitrary function, while u and w are 
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-20 1 ' ' ' ' ' ' ' ' ' 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

y/ 7i 



FIG. 5: 7 as a function of 7 for Examples 2, 3, 4, in arcsec/(year xlO 7 ). 



related to (x,y,z) and to V(t,x,y, z) differently in each 
subfamily. We have 



(r 2 , V) with spherical symmetry, 



2 def a 1 1 1 2 

= .t +y + z ; 



(z, V 7 ) with plane symmetry; 
(x/y, T^/y) with hyperbolic symmetry. 



(8.3) 



The FLRW limit follows when / = and V = 
R(t)g(x,y,z). 

The conformally flat Stephani solution 0, [33[ has the 
metric given by (|8.1[) , the coordinates are comoving, and 
V(t,x,y,z) is given by 



v = h v + \ k{t) [ {x ~ xo{t))2 + {y ~ yo{t))2 



+ (z-zo(t)f]} 



(8.4) 



where (R, k, xq, yo, zq) are arbitrary functions of t. This 
a generalisation of the whole FLRW class, which results 
when (k,xo,yo, zq) are all constant. In general, (I8.4j) has 
no symmetry. 

In these models, in the most general cases, generic null 
geodesies are not RLPs. Consequently, it is not shear 
that causes the non-RLP property. 4 



4 Contrary to what the authors of Ref. Q claim throughout their 
paper. 



In the general type D shearfree normal models, the 
only RLPs are radial null geodesies in the spherical case 
and their analogues in the other two cases. In the most 
general Stephani spacetime, RLPs do not exist. In the 
axially symmetric subcase of the Stephani solution the 
RLPs are those geodesies that intersect the axis of sym- 
metry in every space of constant time. In the spherically-, 
plane- and hyperbolically symmetric subcases, the RLPs 
are the radial geodesies. 

The completely drift-free subcases are conformally flat, 
but more general than FLRW. Their defining property is 
that their time-dependence in the comoving coordinates 
can be factored out, and the cofactor metric is static. The 
FLRW models have the same property. For example, in 
the drift-free spherically symmetric type D case: 



ds 2 = ± { [( A i + A ^ 2 ) PS* dt)] 2 - dr 2 

- r 2 (dtf 2 +r 2 sin 2 tfd^ 2 )} , (8.5) 

the whole non-staticity is contained in V: 

V = B 1 +B 2 r 2 + (A 1 +A 2 r 2 )S(t). (8.6) 

The {Ax, A 2 , B\, B 2 ) are arbitrary constants and S(t) is 
an arbitrary function. This model is more general than 
FLRW because the pressure in it is spatially inhomo- 
geneous. The FLRW limit follows when Ai 7^ and 
B 2 = (A 2 /A 1 )B 1 . 
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IX. DEPENDENCE OF RLPS ON THE 
OBSERVER CONGRUENCE 

The RLPs are defined relative to the congruence of 
worldlines of the observers and light sources. So far, we 
have considered observers and light sources attached to 
the particles of the cosmic medium, whose velocity field 
is defined by the spacetime geometry via the Einstein 
equations. But we could as well consider other timelike 
congruences, or spacetimes in which no preferred time- 
like congruence exists, for example Minkowski. It turns 
out that even in the Minkowski spacetime one can de- 
vise a timelike congruence that will display the non-RLP 
property [34| . 

Take the Minkowski metric in the spherical coordinates 

ds 2 = dt' 2 - dr' 2 - r' 2 (di? 2 + sin 2 M<p 2 ) , (9.1) 

and carry out the following transformation on it: 

t' = (r - t) 2 + l/(r + i) 2 , r' = (r- tf - l/(r + tf. 

(9.2) 



The result is the metric 

ds2= (7T7F {l6u(d<2 ~ dr2) (9 - 3) 

- (u 2 - 1) 2 (dd 2 + sin 2 tfd^ 2 ) } , u d ^ r 2 -t 2 . 

Now we assume that the curves with the unit tangent 
vector field u a = Ur + t) 2 / S a o are world lines of 

test observers and test light sources. 

Proceeding as before we conclude that, with respect to 
this congruence, generic null geodesies in the Minkowski 
spacetime have non-repeatable paths. (The exception are 
those rays that are radial in the coordinates of (|9.4[) ). 
This is because the time-dependence of (|9.4j) cannot be 
factored out. 
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